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A TECHNIQUE FOR SAMPLING POPULATIONS 
OF THE MOUNTAIN PINE BEETLE 


Robert W. Carlson and Walter E. Cole 


Suitable sampling techniques are basic to the determination of relations between an insect 
and the host plant. Development of suitable sampling techniques is part of intensive studies on 
the population dynamics of the mountain pine beetle, Dendroctonus monticolae Hopk., in lodge- 
pole pine, now actively underway in the intermountain area of Utah and Wyoming. 


This paper reports results of research designed to develop a procedure for sampling pop- 
ulations of the bark beetle during the mature larval and pupal stages. The experimental design 
tested for variation between sample sizes and shapes, locations on the tree, and trunk diame- 
ters, as well as testing for the nature of the statistical distribution. 


The study areas selected, though widely separated geographically, were bioclimatically 
equivalent as defined by Hopkins' Law. One was on the Wasatch National Forest in Utah and the 
other on the Teton National Forest in Wyoming. The populations sampled in both areas can be 
considered as epidemic. The Wasatch population has been epidemic in numbers since 1938. 
The Teton population has shown cyclic tendencies of recurrent epidemics; the present cycle 
began about 1958. The affected timber stands differ. The Wasatch stand, in general, consists 
of smaller diameter trees, as old as those in the Teton stand, or older, but less vigorous. The 
Wasatch population was characterized by low brood density and a current static trend. The 
Teton infestation showed a definitely increasing trend. 


EXPERIMENTAL DESIGN 


The experimental work in the two stands was replicated as closely as possible. Three 
variables were measured (counts): (1) density of attack, (2) length of egg gallery, and (3) den- 
sity of brood. All existing stages of brood were recorded: mature larvae, pupae, and callow 
adults. Only successful attacks and consequently successfully constructed egg galleries were 
recorded. 


Sample size: Six sample units were superimposed in a nested fashion at each point sam- 
pled (fig. 1). The six units were: four rectangular--1/10, 1/4, and 1/2 sq. ft.--and a sample 
proportional to tree diameter; and two circular--1/10 and 1/4 sq. ft. Brood counts were made 
only for the rectangular samples. 


Sample location: Each tree was sampled at six points: on the north and south aspects at 
breast height, 5 feet below the height of infestation, and midway between these two points. 


Tree size: Ten trees in each of the following d.b.h. categories per plot were used: 6-9, 
9-12, 12-15, and greater than 15 inches. Only 13 trees per plot were sampled for brood, 
because of the stage of brood development and imminent adult emergence. These trees were 
selected randomly. 
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Figure 1.--The relative placement of the various sample units, 
showing the method of nesting. Scale: 6" =1'. 


ANALYSIS OF DATA 


Statistical distribution: Comparisons of infestations, populations, and other categories 
can be made on the basis of their statistical distributions. In this study, analysis of distribu- 
tions was based upon scatter diagrams of variance plotted over means of the three variables 


measured. 


The scatter of variance in relation to mean attack density and gallery length appeared to 
be random. Comparison of Wasatch and Teton data (figs. 2, 3, 4, and 5) showed no difference 


in distribution. 
ance with an increase in the mean (figs. 6, 7). 


specific contagious distribution involved. 
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The brood data, however, typified a contagious distribution--increasing vari- 
There were not enough data to determine the 


Figure 2. ——letom plot 


--attack density: x 
for) 1/10) (sqz" it. 
rectangular sam- 
ples, o for 1/2 sq. 
ft. samples. 


Figure 3. -- Wasatch 


plot -- attack den- 
sity: x for 1/10 sq. 
ft. rectangular 
samples, o for 1/2 
sq. ft. samples. 


Variance analysis of variance: The assumption of normality must also carry the assump- 
tion that the variances of the means of the variables are homogeneous before parametric meth- 
ods can be used. Like means with like variances are considered to be of one population; 
like means with unlike variances are from different populations. 


The data on attack density and brood were converted to a common basis--1 sq. ft.; 
gallery lengths were converted to a 1/10 sq. ft. basis. 


Complete cross-comparisons of variances were made for attack density and gallery 
length. The 10 trees in each d.b.h. category were considered as replicates (table 1). The 
cross-comparison of variance for brood density was incomplete because the between-plot dif- 
ference was not tested. The difference is obvious because the plots were selected on this - 
basis. The basis for replication for brood density was the 13 trees per plot. 
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Chi-square analyses: Chi-square 
comparisons of variance arrangements 
were used to show the within-tree den- 
sity trends for d.b.h. categories and 
plots. These analyses supplemented 
the variance analysis of variance as a 
tool in making decisions about proper 
sampling procedures. Only data for 
the 1/2 sq. ft. samples were used in 
the analyses. 


Regression analysis: Regression 
and variance analyses were used to 
check the relation of each of the follow- 
ing variables to d.b.h.: height of in- 
festation, attack density, and gallery 


length. 


Figure 6.--Teton plot--brood density: x 


for 1/10 sq. ft. rectangular samples, 
o for 1/2 sq. ft. samples. 


Figure 7.--Wasatch plot--brood density: 
x for 1/10 sq. ft. rectangular samples, 
o for 1/2 sq. ft. samples. 


Table 1.--Form for computation of sums of squares in variance analysis 
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RESULTS 


Sample size: The variance analysis of variances showed significant differences between 
sample sizes, diameters, and plots for the attack density and gallery length variables, but 
showed no significant differences between sample sizes for brood. The Student-Newman-Keul's 
Multiple Range Test (Steel and Torrie 1960)? showed that all sample sizes differed significantly 
from each other for gallery length, and the 1/10 sq. ft. sample differed significantly from all 
others for attack density. The variance between diameter for both attack density and gallery 
length did not show like patterns of occurrence. The test failed to reveal differences between 
sample heights for attack density, and only top and mid-sample heights differed significantly for 
brood density. 


The brood density variation was greater than either the attack density or gallery length 
variation. Thus, inferences drawn from the brood density variable will be used to determine 
sample size for the sampling plan. 


The Teton data showed downward trend of the coefficient of variation with increase of 
sample size, but the Wasatch data did not. This distribution of the coefficient of variation is 
probably due to the effect of zero counts--the larger the sample size the fewer zero counts, and 
the greater likelihood of reducing variance. 


Sample location: The variances of all three variables were generally greater at breast 
height and mid-height than at the top sampling height. The only significant difference was be- 
tween mid- and top heights for the brood density variable. Ona percentile basis, the three 
variables were also more consistent. Variance tends to be greatest at breast height, but since 
the mean values also tend to be larger at this position, the coefficients of variation are not 
correspondingly high. 


Aspect: A significant interaction between aspects by plots indicated that there was never 
complete consistency. In some interactions, as with aspects, the data cannot be treated as if 
they fit the same distribution. Random placement of samples, by aspect, can eliminate this 
elfect. 


Tables 2, 3, and 4 show the results of the variance analysis of variances. 


Tree diameter: Not enough trees were sampled to permit any analysis of brood density 
in relation to d.b.h. 


Regression analysis showed that both attack density and gallery length were significantly 
related to tree size. Attack density and gallery length increased with tree diameter (tables 5, 
6, and 7). 


Size of sample needed: The number of samples (trees) needed for a 20-percent SME at 
the 2/3 probability level was computed for each sample size and all three variables at d.b.h. 
The north and south samples were combined because they were not random with respect to each 
other and in effect constituted a single sample. 


1 Steel, G. D., andj. H. Torrie. Principles and procedures of statistics. New York: 
McGraw-Hill Book Co., Inc. 481 pp. 1960. 


Table 2.--Variance analysis of variances for attack density per square foot + 


Source SUL OE er eeeces oF Mean square EF 
squares : freedom : 

Plots 10668.95 1 10668.95 26.3965" 
Tree sizes 12887 .36 3 4295.79 10..70** 
Sample sizes 107007 .54 5) 21401.51 535.295 
EF 1776.20 3 592.07 1.47 
PS 2491.72 5 498.34 1,24 
TS 7241.00 iS 482.73 1520 
Error. PES 6020.75 15 401.38 
Subtotal 148093 .52 47 
Aspects 0.49 1 0.49 0.00 
Levels 2264.04 2 1132.02 3.267 
AL 411.86 2 205.93 0259 
AP 1530.43 1 1530.43 4,41* 
AT 2087.41 3 695.80 2.00 
AS 1066.13 5 213.23 0.61 
LP SI2Z7 s52 Zz 1863.76 5207" 
ipa 2949.48 6 491.58 1.42 
LS 2327 .64 10 232.76 0.67 
Error, PISAL 72296.93 208 347.58 
Total 236755.45 287 


1Variances computed on a 1 square foot basis. 


Table 3.--Variance analysis of variances for gallery density per 1/10 square foot* 
: Sum of : Degrees of 


S : 
eee squares : freedom Mean eauare : E 
Plots 180.8485 1 180.8485 4.86* 
Tree sizes 8774. 8823 3 2924.9608 78.087" 
Sample sizes 5750.9965 5 1150.1993 30.94** 
Pr 4140.8553 3 1380.2851 sialon 
PS 84.0797 5 16.8159 0.45 
TS 423. 1305 15 28.2087 0.76 
Error, PTS 556.8869 15 37.1258 
Subtotal 19911.6797 47 
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Table 3.--(con.) 


Source Su Degree ok Mean square F 
squares freedom : 

Aspects 1.4056 J) 1.4056 0.00 
Levels 1800.3353 2 900.1676 Woaat 
AL 243.1305 72 12 OOS 0.23 
AP 2786.9334 1 2786 .9334 2.30* 
AT 194.6111 3 64.8704 0.12 
AS 40.5396 5 8.1079 0.02 
LP 119.1520 Ps 59.5760 Oot 
iO 4236.7795 6 706.1299 1.34 
LS 668.6304 10 66.8630 0.13 
Error, PESAL 109753 .6130 208 527.6616 
Total 139756.7831 287 


1 Variances computed on a 1/10 square foot basis. 


Table 4.--Variance analysis of variances for brood density per square foot 


Source 
Plots (Blocks) 


Sample sizes 
Error, PS 


Subtotal 


Aspects 
Levels 

AP 

AS 

AL 

Ihe 

LS 

Error, PSAL 


Total 


Sum of 
squares 


370909825. 
12845676. 
11038342, 


394793845. 


17908167. 
35517526. 
44466345. 

1774537. 
33448518. 
19075248. 

1364835. 
77471934. 


625820956. 


Degrees of 
freedom 
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Mean square 


4281892. 
3679447, 


17908167. 
17758763. 
44466345. 
Sg S12;- 
16724259. 
9537624. 
227472. 
363344, 


Deca 
yall 
1352055 
0.18 
4.96* 
2.83 
0.07 


Table 5.--Regression and variance analysis 
for height of infestation X d.b.h. 


Regression matrix 


vi 2DeE. 
2052.00 
26485.66 
58716.50 80 
1500.00 
6082.70 79 
3685.72 78 


Analysis of variance 


NESS x 
N 80 974.10 
X 12799.57 
xy 
x 938.68 
Yi 
Yi) 
Source: SS 
Bo 52633.80 
By 2396.98 
Res, 3685.72 
Total 58716.50 
x. = D.b.h. 
Y= 


Height of infestation 


MS S F 


2396.98 49,42** 
48.50 


Table 6.--Regression and variance analysis 


for attack density X d.b.h. 
WASATCH PLOT 


Regression matrix 


Nie: X 3 ne 
N 40 487 .50 398.00 
X 6414.29 5160.90 
Y 4454.00 
X 472.88 310.28 
Nd 493.90 
oY 290.31 
Analysis of variance 
sources Di Es: SS : MS 
Bo 1 3960.10 
Br 1 203.59 203.59 
Res. 38 290.31 7.64 
Total 4454.00 
TETON PLOT 
Regression matrix 
IND: X : ¥ 
N 40 486.60 530.00 
X 6385.28 6792.40 
Ms 8016.00 
X 465.79 344.95 
ye 993.50 
ue 738.04 
Analysis of variance 
sources Duk: SS =: MS 
Bo 1 7022.50 
By i 255.46 255.46 
Res. 38 738.04 19.42 
Total 8016.00 
x =)D.b. i. 
Y = Density of attack 
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Table 7.--Regression and variance analysis 
for gallery density X d.b.h. 


WASATCH PLOT 


Regression matrix 


Ni Xx . NC ID) 4 1a 
N 40 487.50 2358.00 
Xx 6414.29 31824.60 
BY 4 169952.00 40 
X 472.88 3086.48 
MG 30947 .90 39 
a: 10802.51 38 
Analysis of variance 
Source: D.F.: SS g MS g F 
Bo 1 139004.10 
Bi 1 20145.39 20145.39 70.86** 
Res. 38 10802.51 284.28 
Total 169952.00 
ELON PLO 
Regression matrix 
IND ee X : Ne SIDA I 
N 40 486.60 PALF 500) 
Xx 6385.28 34359.00 
ny! 203475.00 40 
Xx 465.79 1306.70 
ay 3 18922.78 39 
NE 15257.04 38 
Analysis of variance 
source: De Fs SS 7 MS : F 
Bo 1 184552.22 
Bi 1 3665.74 3665.74 ORT3e* 
Res. 38 15257.04 401.50 
Total 203475 .00 
».G = IDeloialaig 
Y = Density of gallery 


More samples are required for measuring attack density and gallery lengths for the 
smaller size sample than for the larger size. For measurements of brood density, the required 
number of samples decreases as sample size increases. Combining north and south samples at 
d.b.h. generally effected an increase of number of trees needed (table 8). 


Table 8.--The number of trees required to be sampled for a 20-percent SME at 2/3 probability 
level based upon summed north and south bottom samples (rectangular samples only) 


Sample size 


DeneHy ; ait = I /iOlsG tt. 2) 1/4 sae te. “s 1/2 sq.ft. : Proportional 
Attack Teton 9.13 3.36 2.42 ded 
density Wasatch TaI® 4.22 3.63 3.08 
Gallery Teton 6.40 Ser Alt 00 4.67 
density Wasatch 2.46 2.63 2.20 ZeAZ 
Brood Teton 8.19 9.93 8.16 7.56 
density Wasatch 54.06 66.94 67.84 ep el®) 


THE SAMPLING PLAN 


The principal object of any sampling procedure is to secure a sample which, subject to 
physical limitations, will reproduce the characteristics of the population as closely as possible. 
The simplest and most certain way to achieve this goal is for the sample to be drawn at random. 
However, selection at random, subject to restrictions that do not introduce bias--called 
systematic random sampling--is permissible. 


Sampling can be limited to the breast-height region within an infested tree and satisfy the 
requirements of efficiency and reliability in reproducing the characteristics of the population. 
In adhering to the principles of systematic random sampling, the sampling universe was defined 
as that subpopulation within a zone 1 foot above and 1 foot below breast height of infested trees. 


This zone was divided into four quadrants: north, east, south, and west. Ifthe 1/10 sq. 
ft. sample is used, the zone can be divided into six levels of 4 inches each--three above and 
three below breast height; this produces 24 sample loci. If the 1/4 sq. ft. sample is used, the 
zone can be divided into four levels of 6 inches each, producing 16 loci. Prior to sampling, 
two loci were selected at random and sampled without replacement. 


The important advantage of this plan is that an estimate of within-tree variance can be 
computed to determine the statistical distribution; the within-tree error can be treated as error 
attributable to replication; additional degrees of freedom result--one an attribute of replication 
error and the other distributed to the error term or interactions, thus increasing experimental 
precision. This particular two-sample technique will be used in the future work on the popula- 
tion dynamics of the mountain pine beetle in lodgepole pine to determine the type of distribution, 
and to provide the basis for conversion to normality, when necessary. This plan should be a 
step toward the development of more effective sampling techniques. 
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Headquarters for the Intermountain Forest and Range 
Experiment Station are in Ogden, Utah. Project 
headquarters are also at: 


Boise, Idaho 

Bozeman, Montana (in cooperation with Montana 
State University) 

Logan, Utah (in cooperation with Utah State 
University) 

Missoula, Montana (in cooperation with University 
of Montana) 

Moscow, Idaho (in cooperation with the University 
of Idaho) 

Provo, Utah (in cooperation with Brigham Young 
University) 
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